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Abstract 

In this paper we show the existence and form uniqueness of a solution for multidimensional backward stochastic 
differential equations driven by a multidimensional Levy process with moments of all orders. The results are important 
from a pure mathematical point of view as well as in the world of finance: an application to Clark-Ocone and Feynman- 
Kac formulas for multidimensional Levy processes is presented. Moreover, the Feynman-Kac formula and the related 
partial differential integral equations provide an analogue of the famous Black-Scholes partial differential equation 
and thus can be used for the purpose of option pricing in a multidimensional Levy market. 
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1. Introduction 

A linear version of Backward Stochastic Differential Equations (BSDEs in short) driven by Brownian motion 
was initially consided by Bismut(1973)in the context of optimal control. Nonlinear BSDEs were later introduced 
by Pardoux and Peng (1990) and independently by Duffie and Epstein(1992). The BSDE theory has found wide 
applications in partial differential equation theory, stochastic controls and, particularly, mathematical finance(see El 
Karoui and Quenez(1997); Ma and Yong(1999)). 

Situ (1997) studied BSDEs driven by a Brownian motion and a Poisson point process. Ouknine (1998) considered 
BSDEs driven by a Poisson random measure. Nualart and Schoutens (2000) proved a martingale representation 
theorem for Levy processes satisfying some exponential moment condition. By using this martingale representation 
result, Nualart and Schoutens (2001) established the existence and uniqueness of solutions for BSDEs driven by a 
Levy process of the kind considered in Nualart and Schoutens (2000), furthermore Nualart and Schoutens (2001) 
presented the Clark-Ocone and the Feynman-Kac formulas, the related Partial Differential Integral Equation (PDIE) 
and their applications in finance. 

In the past twenty years, there is already a growing interest for multidimensional Levy Processes. Some con- 
cepts and basic properties about multidimensional Levy Processes were summarized in Sato (1999). Applications 
of multidimensional Levy Processes to analyzing biomolecular (DNA and protein) data and one-server light traffic 
queues were explored by Dembo, Karlin and Zeittouni (1994). A small deviations property of multidimensional Levy 
Processes were discussed by Simon (2003). In finance research, practically all financial applications require a multidi- 
mensional model with dependence between components: examples are basket option pricing, portfolio optimization, 
simulation of risk scenarios for portfolios. In most of these applications, jumps in the price process must be taken into 
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account. Cont and Tankov (2004) systematically investigated these problems in multidimensional Levy market. In 
addition, the optimal portfolios in multidimensional Levy market is discussed by Emmer and Kliippelberg (2004), and 
option pricing is investigated by Reich,N., Schwab, C. and Winter, C.(2009). Some simulation approaches for mul- 
tidimensional Levy processes are also investigated in Cohen and Rosihski (2007). Levy copulas was also suggested 
by Kallsen and Tankov (2006) in order to characterize the dependence among components of multidimensional Levy 
Processes. 

Recently, a martingale representation theorem for multidimensional Levy processes was also proved in Lin (201 1), 
and the obtained representation formula was similar as that in Nualart and Schoutens (2000). The purpose of this paper 
is to use this martingale representation result obtained in Lin (201 1) to establish the existence and form uniqueness of 
solutions for BSDE's driven by a multivariate Levy process considered in Lin (201 1). Although the proof techniques 
are similar to those in Nualart and Schoutens (2001), the results are important from a pure mathematical point of view 
as well as in the world of finance. This is illustrated in the applications. The resulting Clark-Ocone and Feynman-Kac 
formulas are fundamental ingredients in the build up of an Malliavin calculus for multidimensional Levy processes. 
Moreover, the Feynman-Kac formula and the related Partial Differential Integral Equation (PDIE) also have an im- 
portant application in finance: they provide us an analogue of the famous Black-Scholes partial differential equation 
and is used for the purpose of option pricing in a multidimensional Levy market. 

The paper is organized as follows. Section 2 contains some preliminaries on multidimensional Levy processes. 
Section 3 contains the main result on BSDEs driven by multidimensional Levy processes. In Section 4 we have 
included some applications of BSDE's driven by multidimensional Levy processes to the Clark-Ocone, the Feynman- 
Kac formulas, and option pricing in a multivariate Levy market. Finally, in the appendix one can find detailed proofs 
of the main results. 

2. Preliminary 

A Revalued stochastic process X = [X(t) = (Xi(t),X 2 (t), ■ ■ ■ ,X n (t))',t > 0} defined in complete probability 
space (Q, P) is called Levy process if X has stationary and independent increments and X(0) = 0. A Levy 
process possesses a cadlag modification and we will always assume that we are using this cadlag version. If we let 
& t = %M ,jV, where % = cr{X(s), < s < t) is the natural filtration of X, and jY are the P-null sets of J? then 
{^r, t > 0} is a right continuous family of cr-fields. We assume that & is generated by X. For an up-to-date and 
comprehensive account of Levy processes we refer the reader to Bertoin (1996) and Sato (1999). 

Let X be a Levy process and denote by 

X(t-) = lim X(s), t > 0, 

s— >t,s<t 

the left limit process and by AX(t) = X(f) - X(t-) the jump size at time t . It is known that the law of X(f) is infinitely 
divisible with characteristic function of the form 

E [exp(i0 ■ X(t)j] = (0(0))' , = (6>i, 2 , • ■ • , e n ) e W 

where <p(0) is the characteristic function of X(l). The function tf/(0) = log<p(0) is called the characteristic exponent 
and it satisfies the following famous Levy-Khintchine formula (Bertoin, 1996): 

t//(0) = -l-0-Ifi + ia-0+ \ (exp(i0-x)-\-i0-x\\ x \<i)v(dx). 

where a, x e W, 2 is a symmetric nonnegative-definite n x n matrix, and v is a measure on M"\{o} with J (\\x\\ 2 A 
l)v(dx) < oo. The measure v is called the Levy measure of X. 

Throughout this paper, we will use the standard multi-index notation. We denote by No the set of nonnegative 
integers. A multi-index is usually denoted by p, p = (p\,pi, • ■ ■ ,p n ) e Nj>. Whenever p appears with subscript or 
superscript, it means a multi-index. In this spirit, for example, for x = (x\ , ■ ■ ■ , x n ), a monomial in variables x\ , • • ■ , x n 
is denoted by x p = x 1 ^ ■ ■ ■ x%" . In addition, we also define p\ = p\ ! • • ■ p n \ and \p\ = p 1 + ■ ■ ■ + p n \ and if p, q e N" y 
then we define S p ^ q = 6 Puqi ■ ■ ■ (5 p „, g „. 
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Hypothesis 1. We will suppose in the remaining of the paper that the Levy measure satisfies for some e > 0, and 
A>0, 



f ■ 

J\x\>e 



exp(A\\x\\)v(dx) < oo. 



This implies that 



x p v(dx) <cv. \p\>2 



and that the characteristic function E [exp(iO ■ X(t))] is analytic in a neighborhood of origin o. As a consequence, X(t) 
has moments of all orders and the polynomials are dense in L 2 (K", P o X(f) -1 ) for all t > 0. 

Fix a time interval [0, T] and set L\ = L 2 (£l, & T , P). We will denote by & the predictable sub-cr-field of & T ® 
&[0,t]- First we introduce some notation: 

• : Let Hj denote the space of square integrable and ^-progressively one-dimensional measurable processes 
4> = {(f>(t),te [0, T]} such that 



-LT 



U(t)\\ 2 dt 



< CO. 



• : Ml. will denote the subspace of Hi formed by predictable processes. 



• : (Hj(l )) m and (Mj(l )) m are the corresponding spaces of m-dimensional r-valued processes equipped with 
the norm 



E 



Jr-T 00 
d=\ peW, 



1,2, • • • ,m, 



\w\ 



= Yj\\Mt)\\ 2 P , 



k=l 



def 



where = (4> u <t>2, • ■ ■ , 0m)', «t = : p e Nq}, it = 1,2, ■ • ■ , m and N» = {p e : |p| = d}. 
• : Set ^=ffjx {M 2 T {l 1 )) m . 

Following Lin (201 1) we introduce power jump monomial processes of the form 



X(t) 



(Pl.-.Pn) 



def 



0<s<t 



The number \p\ is called the total degree of X(t) p . Furthermore define 



Y(t) 



(pi,-,p„) 



def 



X(t) 



(Pl,~,Pn) 



E[X(t) 



(Pl,-,Pn) 



] = x(ty 



(P\,-,Pn) 



- m„t, 



the compensated powerjump process of multi-index/? = (p\,p2,--- , p„)wheiem p = J ]~[ x f v (^)- Under hypothesis 

i=i 

1, Y(t) {pu '"' p, < ) is a normal martingale, since for an integrable Levy process Z, the process [Z, - E[Z,],t > 0} is a 
martingale. We call Y(fy Pl, '" ,p '" the Teugels martingale monomial of multi-index (p\, ■ ■ ■ ,p n ). 

We can apply the standard Gram-Schmidt process with the graded lexicographical order to generate a biorthogonal 
basis {H p ,p e N"}, such that each H p (\p\ = d) is a linear combination of the Y q , with |^| < \p\ and the leading 
coefficient equal to 1 . We set 

W-i 

H p = Y p + J] c q Y« + Y,Y, c « Y "> 

q<p,\q\=\p\ k=l \q\=k 
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where p = [p\, ■ ■ ■ ,p n ), q = [q\, ■ ■ ■ ,q n ] and -< represent the relation of graded lexicographical order between two 
multi-indexes. Some details about the technique and theory of orthogonal polynomials of several variables refer to 
Dunkland Xu (2001). 



set 



Set 



P(x) p 
P(x) p 



= x" + 



= x p 



2 

q<p,\q\=\p\ 



CnX" + 



CqX^ , 



q<p,\q\=\p\ 



k=\ \q\=k 
\p\-l 

k=2 \q\=k 



c g x", 



H p (t) 



Specially we have 



H ei (t) 
H e2 (t) 

H e "(t) 



2 (AX 1 r---(AX n y»+ J] c ?( A *i) 9 

0<s<t V q<P,\q\=\p\ 
lp|-l 

+ ^^C 9 (AX 1 ) <?I ---(AX„)^ 



•(AX„) 



*=1 l«l=* 
-tE 



9<P,I?I=W *=1 l«l=* 

(c ei X!(l) + • • • + C e X n {\)) + £ p(AX(i)) 



0<.s<f 



-fE 



X P(aX(s)) 



0<s<; 



-?E[c ei Xi(l) + --- + c e X„(l)]. 



Q.dXX^O-fE^d))), 

c. 2 (2)(X 2 (r) - fE(X 2 (l))) + c ei (2)(X 1 (f) - fE(X!(l))), 

c e »(X„« - WX n (l))) + c,^(«)(^i(0 - tE^n-itf))) 
+ --- + c ei (ri)(Xi(f)-WXi(l))). 



(1) 



The main results in Lin (2011) is the Predictable Representation Property (PRP): Every random variable F in 
L 2 (Q, has a representation of the form 

CO rj. 

F = E(F) +ni <$> p (s)dH p {s) 

where O^s) is predictable. It is worthwhile to emphasize that <& p (s) is not uniqueness, that is to say that <b p (s) is 
different for different Gram-Schmidt process, but the form of & p (s) is uniqueness. This type of uniqueness is called 
"form uniqueness". This result is an extended version for the corresponding Theorem in Nualart and Schouten (2000). 

Remark 1. If v = 0, we are in the classical Brownian case and H p (t) = 0, \p\ > 2. If p has only mass in 1, we are in 
the Poisson case; and also here H p (t) = 0, \p\ > 2. Both case are degenerate cases in this Levy framework. 

From these observations, it is not so hard to see that the PRP property shows that financial markets based on 
a non-Brownian or non-Poissonian Levy process, i.e. with a stock price behaviour S,-(r) = S i(0)exp(rt + X,(f)), 
i = 1, 2, • • • ,n, are so called incomplete, meaning that perfectly replicating or hedging strategies do not exists for all 
relevant contingent claims. 
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3. Multidimensional BSDEs Driven by Multidimensional Levy Processes 

Taking into account the results and notation presented in the previous section, it seems natural to consider the 

CO 

-dY(t) = /[aHZ(p-J^(« Y(T) = g, (2) 

d=l peN'' 

where 

. :Y(f) = (Y 1 (f),Y 2 (f),--- ,YJf)y. 

• : Z(t) = {z p (t)}pew , each component z p (t) = (z p , ■ ■ ■ ,Zm)' is a m-variables & t predictable function; 

• f = (fufi, - ■ ■ , fm)' : Q X [0, 71 x W" x (m 2 (/ 2 )) m -» W n is a measurable m-dimensional vector function such 
that /(-,0,0) E (// 2 ) m 

• f is uniformly Lipschitz in the first two components, i.e., there exists Ct > 0, k — 1, 2, • • • , m, such that dt ® dP 
a.s., for all (y u zi) and (y 2 , Zi) in M ffl x (Z 2 ) m 

\f k (t,yu zi) - f k (t,y 2 , z 2 )\ < C k (||vi - y 2 \h + llzi - z 2 || ( p r ) , k = 1, 2, • • • ,m. 
. fei 2 r (0,P). 

If (/,£) satisfies the above assumptions, the pair (/,£) is said to be standard data for BSDE. A solution of the 
BSDE is a pair of processes, {(Y(t), Z(f)), < t < T) e H\ x [m\ (Z 2 ))'" such that the following relation holds for all 
fE[0,r]: 

Y(f)=€+ f f(s,Y(s-),Z(s))ds-J]J] [ z p (s)dHP(s). (3) 

Note that the progressive measurability of {(Y(t), Z(t)), < t < T) implies that (T(0), Z(0)) is deterministic. 
When the orthogonal polynomials are fixed, A first key-result concerns the existence uniqueness of solution of 
BSDE: 

Theorem 1. Given standard data (/, £), there exists a unique form solution (Y, U, Z) which solves the BSDE (3) 

The proof can be found in the Appendix, as the proof of the continuous dependency of the solution on the final data £ 
and the function /. 

Theorem 2. Given standard data (/, and (/', %'), let (Y, Z) and (Y', Z') be the unique adapted form solutions of 
the BSDE(3) corresponding to (/, and (/', £'). Then 



E 



\ 

ds 



f \\Y(s-) - Y(s-)\\ 2 + £ Z - ~ z "^\\ 2 

Jo { d=lpeN» 

< C (e[||£ - 1|| 2 ] + E [JT \\f(s, Y(s-), Z(s)) - f(s, Y(s-), Z(s))\\ 2 ds J . 

The definition of "form solution" in Theorem 1 and Theorem 2 is deferred until later in the proofs of Theorem 1 
and Theorem 2 in Appendix. 



4. Applications 

Suppose our n-dimensional Levy process X(t) has no Brownian part, i.e. X(t) = at + L(t), where a = (a\, ■ ■ ■ , a n )' 
and L(f) is n-dimensional pure jump process with Levy measure v(dx). 
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4.1. Clark-Ocone Formula and Feynman-Kac Formula 

Let us consider the simple case of a BSDE where / = 0, and the terminal random vector £ is a function of X(T), 
that is, 



• dY(t) = - £ 2 *WH'{t), Y(T) = g(X(T)) 



or equivalently 



F(f) = g(X(T)) - > ' > ' I z p (s)dH?(s), (4) 



d=l neN" ^' 



where g = (gi,g2, ■ ■ ■ >&»)' and E(||g(X(r))|| 2 ) < oo. Let6 k = 6 k (t,x),k = 1,2, • ■ ■ ,m, be the solution of the following 
PDIE(Partial Differential Integral Equation) with terminal value g^: 

^(t,x)+ f ie k (t,x + y)-e k (t,x)-^(t,x)-y)v(dy) + a-^(t,x) = 0, 
at J R „\ dx ) dx 

6 k (T,x) = g k (x), (5) 

where a = a + f yv(dy). Set 



%\t,x,j) = 9 k (t,x+y) - 9 k (t,x) - |£(f, x) y. (6) 



The following result is a version of the Clark-Ocone formula for functions of a Levy process. Again the proof can be 
found in the Appendix. 

Proposition 3. Suppose that 6 k is a C 1 ' 2 function for k — 1 , 2, • • • , m such that || and |Jr ® re bounded by a polyno- 
mial function ofx, uniformly in t, then the unique adapted form solution of (4) is given by 

Y k (t) = e k (t,x(t)) 

4 = f M J ( k ) (t,X(t-),y)pP(y)v(dy), \p\>2, 

= Jk X(t-),y)p e >(y)v(dy) + £ f*c (7 , i = 1,2, • • ■ ,n 

7=1 ' 



where 6 k = 6 k {t,x)for k — 1,2, • • • ,m are f/ze solutions of the system of PDIE(5), (f k \t,x,y) for k — 1,2, • • • ,m are 
g/ven fry (6) ana" [c,j] is f/ze inverse matrix of the coefficient matrix Gram-Schmidt transformation (1). 

Now by taking expectations we derive that the solution 0(t,x) = (6i(t,x), • • • , 6 m (t,x))' to our PDIE(5) equation has 
the stochastic representation 

0(t,x) = E[g(X(T))\X(t)=x]. 

This is an extension of the classical Feynman-Kac Formula. 

If Ln \\y\\ v (dy) < oo, and we take a = L „ ,,yv(dy), then the equation (5) reduces to 



™{t,x)i \ (0U.x+y)-0u.x))v(dy) = 0. 



0(7\x) = g(x), 

and we have for k = 1 , 2, • • • , m 

< = f R M(t,X(t-)+y)-0 k (t,X(t-)))pP(y)v(dy), \p\>2, 

z* = + J) " WU«-MM4f) + I $Cij, i = 1,2, ■ • • ,n. 

7=1 ' 
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Example: In this example, we define a two-dimensional Poisson process by using Levy copulas. All the concepts and 
notations are adopted from the Kallsen and Tankov (2006). Here the two marginal processes are given respectively by 
two Poisson process Ni(t) where i = 1, 2. In particular, the Levy copula F(u\,ui) : R" — > IR is taken as 

F(u u ui) = (\ur\-" + taP) -1 '" (^{«,« 2 >o} - (1 - ^)/{«,« 2 <o|) (7) 

where p > and any n e [0, 1 ] . It defines a two parameter family of Levy copulas which resembles the Clayton family 
of ordinary copulas. Thus its Levy measure can be calculated as: 

v(dxidx 2 ) = ;— 6i(xi)6i(x 2 ). (8) 

Set two compensated Poisson process X,-(f) = Ni(t) - Ait for ;' = 1,2, then we can construct a set of martingales 
H p , p e NjJ, i = 1, 2 by orthogonalizing procedure proposed by Lin (201 1) such that H p is strongly orthogonal to H 9 
with respect to Levy measure v(dx\dxi), for p + q. Moreover the PDIE (5) reduces to 

+ m.x + 1) - (Kt,x)) - %,x) ■ l ^+^f = 0, 

9(T,x) = g(x), 

The Clark-Ocone Formula is now given by 

n „t 

g k (X(T)) = E[g k (X(T))] + J] (0 k (s,X(s) + l)-e k (s,X(s)))dX i (s) 

4.2. Nonlinear Clark-Haussman-Ocone Formula and Feynman-Kac Formula 
Let us consider the BSDE 

CO 

-dY(t) = f(t,Y(t),Z(t))-J] J] Z p {f)dH p {f), Y(T) = g(X(T)) (9) 

d=\ peN" 



or equivalently 



XT *> pT 

f(s,Y(s-),Z(s))ds-J]J] Z p (s)dHP(s), 
d=l peN" J' 



at J R » 



Suppose that 6 k = Q k {t , x) for = 1 , 2, • • • , m satisfy the following system of PDIE: 

fli'V.*.^) + « ' ?n('.*) + /* (t,e k (t,x),G k (t,x)) = 0, (10) 
ax 

6 k (T,x) = gk(x). 

where as in the previous section, we define 9f\t,x,y) by (6), and ® k (t) = {0^(t)} pe N«, where 

tf(t,x) = f 6 { '\t,x,y)f(y)v{dy), (11) 



Proposition 4. Suppose that 6 k is a C 1 ' 2 function such that ^| and are bounded by a polynomial function of x, 
uniformly in t, then the unique adapted form solution of (9) is given by 



Y k (t) = 6 k (t,X(t)) 
z p = f 



- f R J k l) (t,X(t-),y)p»(y)v(dy), \p\>2, 



JEp (f k l \t, X(t-),y)p^(y)v(dy) + £ £c, 7 , i = 1,2, ■ • ■ ,n. 

7=1 ' 



where 9 k = 9 k (t,x)fork = 1,2, ■•■ ,m are the solution of the system of PDIE( 10), 6^\t,x,y) is given by (6) and[dij] 
is the inverse matrix of the coefficient matrix Gram-Schmidt transformation (1). 
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Notice that taking expectations we get 



6 k (t, x) = E [g k (X(T))\X(t) = x] + E [ f /* (s, k (s, X(s-)), ® k (s, X(s-))) ds\X{i) = x . 



t 



Example: Consider again the very special case where have a two-dimensional Poisson process (Ni(t),N2(tJ) by using 
Levy copulas (7). Set X t (t) = N,(f) - A t t where i = 1,2. Then the PDIE (10) reduces to 



where k = 1,2, • • • ,m. 
4.3. Option Pricing 

In the last two decades several particular choices for one-dimensional non-Brownian Levy processes were pro- 
posed. Madan and Seneta (1990) have proposed a Levy process with variance gamma distributed increments. We 
mention also the Hyperbolic Model proposed by Eberlein and Keller (1995). In the same year Barndorff-Nielsen 
(1995) proposed the normal inverse Gaussian Levy process. The CMGY model was also introduced in Carr et al. 
(2000). Finally, we mention the Meixner model (see Grigelionis (1999) and Schoutens (2001)). All models give a 
much better fit to the data and lead to an improvement with respect to the Black-Scholes model. 

Multidimensional models with jumps are more difficult to construct than one-dimensional ones. A simple method 
to introduce jumps into a multidimensional model is to take a multivariate Brownian motion and time change it 
with a univariate subordinator (refer to Cont and Tankov (2004)). The multidimensional versions of the models 
include variance gamma, normal inverse Gaussian and generalized hyperbolic processes. The principal advantage 
of this method is its simplicity and analytic tractability; in particular, processes of this type are easy to simulate. 
Another method to introduce jumps into a multidimensional model is so-called method of Levy copulas proposed by 
Kallsen and Tankov (2006). The principle advantage in this way lies in that the dependence among components of 
the multidimensional Levy processes can be completely characterized with a Levy copula. This allows us to give a 
systematic method to construct multidimensional Levy processes with specified dependence. 

Here we define a multivariate Meixner process by using Levy copulas, and all the concepts and notations are 
adopted from the Kallsen and Tankov (2006). In particular, for n > 2, the Levy copula F{u\, ■ ■ ■ , u„) : R" — » R is 
taken as 



It defines a two parameter family of Levy copulas which resembles the Clayton family of ordinary copulas. It is in 
fact a Levy copula homogeneous of order 1 , for any fi > and any 77 € [0, 1] . 

In addition, we know that if the tail integrals £/,(x,) related to the marginal Levy densities Vj(dxj), ; = 1, • • • , n, are 
absolutely continuous, we can compute the Levy density of the Levy copula process by differentiation as follows: 




e k (T,x) = gk (x), 

k= 1,2,- ■■ ,m 



And we derive the nonlinear Feynman-Kac Formula: 




•T 



f k (s, 0k(s, X(s-)), ® k (s, X(s-))) ds\X(t) = x . 




v(dxi,--- ,dx n ) = d\ ■■■d n F\^ =Ul(xil ...^ n=Un(Xn) vi(dx l )---v i _(dx n ) 
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where v\(dx\), • • • , v n (x n ) are marginal Levy densities. 

Hence we are able to construct a n-variate Levy process X(f) = (X\(t),X2(t), ■ ■ ■ ,X n (t))' with state space W and 
characteristic triple (P, vp, a) by using the above complete dependent Levy copula or Clayton Levy copulas proposed 
by Kallsen and Tankov (2006). Here, in particular, we will consider Levy marginal measure v, of Meixner type for 
Levy copula. We recall the marginal density fi, the cumulative generating function Kj, the drift a„ and the marginal 
Levy measure v,, for the ith component of the Meixner Process (X,(f), t>0},i= 1, 2, • • • , n, are illustrated as follows, 

- fi{Xi,ai,Pi,8i,Hi) = 



KiiOc, an,Pi, 8 h = ufa + 2£, [logcos^ - logcos " 1 ' ) , 

8 x- 

«. r"° sinh^ 

*J 1 rv. 



Vi(dxi;ai,/3i,Si,Hi) = ' 6 dx h 
XiSinn— 

where a t > 0, -n < /?, < n, e K, and 5, > 0. 

From the form of the cumulative generating function one easily deduces that the density at any time t can be 
calculated by multiplying the parameters 6, and by t for both cases. 

In addition, Reich, etc. (2009) had proved that e Xi is a martingale with respect to the canonical filtration & of X if 
only if 



L Jr» 



(e^-l- Zj I lzl < i )v(dz) = 0. 

Following Reich, etc. (2009), we assume a risk-neutral market which consists of one riskless asset (the bond) with 
price process given by B(t) = e rt , where r is compound interest rate, and n(> 1) risky assets (the stocks), with price 
process: 

S ,(f) = S i(0)exp(rt + X^t)), i = 1, 2, • • ■ , n, 

where X(t) = (X\(t), X2(t), • • • , X n (t))' is a n-variate Levy process and characteristic triple (P, vp, a) under a risk- 
neutral measure P such that e x > is a martingale with respect to the canonical filtration &® = f <t(X(s), s < t), t > 0, of 
the multivariate process X. Denote by F(dx) the probability measure of X(l). 

We consider a European option with maturity T < oo and payoff G(S) where 5 = (S i, • • • ,S„Y, and we assume 
that G(S) is Lipschitz. According to the fundamental theorem of asset pricing (see Delbaen and Schachermayer 
(1994)) the value V(t, S) of this option is given by 

V(t,S) = E[e- r(:r - G(S(r))|S(f) = S]. 

If V(t,S) satisfies 

V(t, S) e C 1 ' 2 ((0, T) x R n >Q ) n C° ([0, T] x R£ ) 
Then Reich,etc.(2009) had proved that V(t, S) is a classical solution of the backward Kolmogorov equation: 



f ( t , S) + \ p^o-^i.S) + rp t %.«,S) - rVit,S) 
+ jT \V(t,S ■ e z ) - V(t,S) " Z S '<«* ~ l )^-(t,S)Jv(dz) = 0, 

def 

on (0, T) x M" where V(t, S ■ e z ) = V(t, S \e Zl , • • • , S n e z "), and the terminal condition is given by 

V(T,S) = g(S), VSeR" £0 . 
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5. Appendix: Proofs of the Results 



Proof of Theorem 1: 

We define a mapping O from 'W 2 into itself such that (Y, Z) e Hj is a solution of the BSDE if only if it is a fixed 
point of O. Given (U, V) e ( H\, we define (Y, Z) = 0>([7, V) as follows: 



F(0 = £ f+ /(j,l/(j-),V(j))dj|f; 



< / < T, 



and (Z(f), < t < T] is given by the martingale representation of Lin (201 1) applied to the square integrable random 
variable 

f /(*,£/(*-), F(*))J S , 
Jo 

i.e., 

" - £ f(s,U(s-),V(s))ds 

£ + £f(s,U(s-),V(s))ds]+Z 2 JT r cf>P(s)dHP(s), 

In according to the result in Lin (2011), a different Gram-Schmidt process will generate a different <f> p (s). Here 
and hereafter, for a given (U, V) e Hj, we always mean that the <& p (s) is generated by arbitrarily selecting a fixed 
Gram-Schmidt process. It is to say that the Gram-Schmidt process is fixed once the Gram-Schmidt process is first 
arbitrarily selected. In the following, the existence and uniqueness of the solution for BSDE(3) will related to this 
selected W(s). Although a <£> p (s) is different for a different Gram-Schmidt process, the expression forms of W(s) are 
the same, the corresponding solutions are called "form solution". 

Taking the conditional expectation with respect to f, in the last identity yields 



Y, + f(s, U(s-), V(s))ds = To + X X 4>»( S )dHP(s), 



from which we deduce that 

XT °° pT 

f(s,U(s-),V(s))ds-J]J] <pP(s)dHP(s), 
d=\ P eW d J < 

and we have shown that (Y, Z) e 1-Oj. solves our BSDE if only if it is a fixed point of (J). 
Next we prove that <5 is a strict contraction on 'Hj equipped with the norm 

aT vl/2 
e*(||r(s-)|| 2 + ||Z(s)|| 2 )£fol , 

for a suitable /? > 0. Let (U, V) and (£/', V) be two elements ofHj and set <D(t/, V) = (Y, Z) and <D(t/', V) = (F, Z'). 
Denote (T7, V) = (1/ - V , V - V) and (Y,Z) = (Y - Y' , Z - Z'). 

Applying Ito's formula from s = t to s = T, to e^ J ||T(i) - T(i)'|| 2 , it follows that 

^'||F(f) - T(f)'ll 2 = - F'(i-)ll 2 di 

-2 J - • d(Y(s) - Y(s)') 

"ZJ ^ S d[Yj-Y' p Yj-Y' } ](s). (12) 
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We have 



■d(Y(t)-Y(tY) = (f(t,U(t-),V(t))-f(t,U(t-)',V(ty))dt 

CO 



dVYj-Y'yYj-Y'fif) 

= hh^#w ff '* 7 = 1,2, ■•■,»! 

where the symbol "o" represents the Hadamard-Schur product for two vectors. 

<H",H 9 >(t) = 6 M t. 

Hence, taking expectations in (12), we have 



E \j*\\Y(f) - Y(t)'\\ 2 ] +ZZ P [ ef'tfisfds 
L J j=ld=lpeW d v jt 

f^\\Y(s-)-Y'(s-)\\ 2 ds 



-ySE 

+2E 



/ - Y'(s-)) ■ (f(s, U(s-), V(s)) - f(s, U(s-)', V(s)'))ds 



Using the fact that / is Lipschitz with constant C yields 

rT „„.,— 



E[e^||y(0-F(f)'ll 2 ]+E JT e^||0(i)|| 2 ^ 

+2CE[j;V s ||y( i -) - F'(5-)|| ■ (\U(s-) - U'(s-)\ + \\V(s)\\)ds 

If we now use the fact that for every c > and a, b e R we have that 2ab < ca 2 + -b 2 and (a + b) 2 < 2a 2 + 2b 2 , we 
obtain 

" rT fl„„-r. 



E [e*||r(0 - F(f)'ll 2 ] + E [JT e^\\4>(s)\\ 2 ds 
fef"\\Y(s-)-Y'(s-)\\ 2 ds 



< (4C 2 -/?)E 

+ ^(\U( S -) - U'{s-)\ 2 + W(s)\\ 2 )ds 

Taking now /3 = AC 2 + 1, and noting that eP*E [\\Y(t) - F(f)'ll 2 ] > 0, we finally derive 

E J t T eP s \\Y{s) ~ Y(s)'\\ 2 ds\ + E [j^ eP s \^{s)\\ 2 ds 
\ f ^ (\U(s-) - U'(s-)\ 2 + || V(s)|| 2 ) ds 



\E 



that is 



ll(F,Z)|| 2 < -\\(U,V)\\j, 

from which it follows that <5 is a strict contraction on 'H 2 , equipped with the norm || • 1^ if = AC 2 + 1. Then O has a 
unique fixed point and the theorem is proved .o 
Proof of Theorem 2: 
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Applying Ito's formula from s = t to s = T, to \\Y(s) - F(s)'|r, it follows that 

\\Y(T) - Y(T)'\\ 2 - \\Y(t) - Y(t)'\\ 2 = 2 J - Y'(s-)) ■ d(Y(s) - Y(s)') 

m „T 

+ 2 d[Yj-Y>,Yj-Y>](s). 

Taking expectations and using the relations 

-d(Y{t)-Y{t)') = {f{t,Y(t-),Z{t))-f'(t,Y(t-)',Z(t)'))dt 

oo 

+YjTj* P(t)dH " {t) ' 

d=\ peN) 



diYj-Y'^Yj-Y'^t) 



CO oo 



2 2 2 2 /j(t)4>"j(t)d[H p ,H^(t), j=\,2,---,m 



we have 



<H",H q > (f) = <W- 



l|y(0 - y(0'll 2 + 2 2 2 e X *'(*) 2 <fr 



= E 



:[X r i^-ni 

+2E [ J t T (Y(s-) - Y'(s-)) ■ (f(s, Y(s-), Z(s)) - f'(s, Y(s-)', Z{s)'))ds 
Using the Lipschitz property of /', and computations similar to those of the proof of Theorem 1 we obtain 

Ef||F(f)-y(f)'ll 2 ] + 5E[j; r ||0WII 2 ^ 
< E [|£ - ? || 2 ] + (1 + 2C + 2C' 2 )E [j; 7 \\Y(s-) - Y'(s-)\\ 2 ds 
+E \\f(s, Y(s-), Z(j)) - /'(j, r(j-)', Z(s)')|| 2 ^v] . 

Then by Gronwall's inequality the result follows. o 

Lemma 5. Lef /i:flx[0, r]xR"->Rka random function measurable with respect toV® such that 

\h(s,y)\<a s (yyA\\y\\) a.s., (13) 

where [a s , < s < T] is a nonnegative predictable process such that E a 2 s ds < oo. Then for each t e [0, T] we 
have 

XT f °° r-t 

h{s,y)v{dy)ds + Y J Yu <K.s,-),f > dH>{s). 

Proof of Lemma 5: Because (13) implies that E [jTf J Rrl \h(s, y)\ 2 v(dy)ds] < oo, we have that 

M(t) = V h(s,AX(s))- f f h(s,y)v(dy)ds. 
o7T<t Jo Jr " 
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is a square integrable martingale. By the Predictable Representation Theorem, there exists a process <f> in the space 
(M 2 T (l 2 )) m such that 

CO , 

M{t)= 2 2 f cf>P(s)dHP(s) 

d=\ peWj 

Taking into account that < H p , H q >,- tS pq , we have 

<M,H p >,= f <f> p (s)ds. (14) 
Jo 

On the other hand, using that aM(s)aH p (s) = h(s, AX(s))pP(AX(s)) we obtain 

<M,H">,= f f H.s,y)?{y)v(dy)ds. (15) 
Jo Jr» 

Consequently, (14) and (15) imply 

^(j)= f ^.jO^M^). 

and the result follows. o 
Proof of Proposition 3: 

Under the hypotheses of Proposition 3 the function 0^\t,x,y) for k = 1,2, ••• ,m given by (6) satisfies the 
hypotheses in Lemma 5 imposed on h due to the mean value theorem, when we take x = X(t-). 
Apply Ito's lemma to 6 k (s, X{s)) from s = t to s = T: 

e k (T,X(T))-6 k (t,X(t)) = J ^(s,X(s-))ds + J]£ ^-(s^is-VdXiis) (16) 

+ J] Ws,X(s))-6As,X(s-))-^^(s,X(s-))AX,is) . 

t<s<Tl i=l X ' 

n „ 

If we apply Lemma 5 to h(s, X,y) = 6 k (s, X(s-) + v) - 6 k (s, X(s-)) - 2 iris, X(s-))yi, we obtain 



2 

1<S<T 



e k (s,x(s)) - e k (s,x( S -)) - 2 ^xis-VAXjis) 

i=\ 



= 22 f(J R J k 1) (s,X(s-),y)pP(y)v(dy))dHP(s) ( 17 ) 

d=l peNJ 

+ f f Rn ^ k \s,X(s-),y)v(dy)ds 
Hence, substituting (17) into (16) yields 



gk (X(T))-e k (t,x(t)) 

= f f(s, X(s-))ds +Zf %(s, Xis-VdXJs) 

i-l 1 

+ 2 2 f t T (f R „e ( k l) (s,X(s-),y)pO(y)v(dyj)dHP(s) 
+ J t T i n O ( k \s,X(s-),y)y(dy)ds 

Notice that 

n 

Xi(t) = Yf\t) + tE(Xi(i)) = 2 ~ c a Hei W + ^(i)). 

7=1 
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(18) 



and 



E(X i (l)) = a i + f y iV (dy). 
By applying the condition (5) and F(0) = E[F(0)] = E[g(X(T))], We can rewrite (18) as 
g k (X(T)) = E[g k (X(T))] + y f ^(j,jr(j-))Ycydff*>(0 

1=1 J ' J=l 
+ ZZ f (f < 1) (a(H,j)/(y)vw)ffi ? W (19) 

which completes the proof of the Proposition, o 
Proof of Proposition 4: 

Apply Ito's lemma to 9 k (s, X(s)) for k = 1,2, • • • ,d from s = t to s = T. By using Lemma 5, we obtain the 
equality (18). Now, using (10) we get 

g k (X(T)) - k (t,X(t)) = - ff k (s, Y(s-),Z(s))ds 

+ if d £(s,X(s-)) £ CijdH'Ks) + £ f(J R „ ( k l \s,X(s-),y)p^(y)v(dy))dH e >(s) 

i=l ' j=\ i=\ 

+ h J t T (f Rn e^(s,X(s-),y)pP(y)v(dy))dHP(s) 

d=2 peWj 

which completes the proof of the Proposition. 
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